Rashba interaction in quantum wires with in-plane magnetic fields 



in 
o 
o 

(N 

o 

Q 



Lloreng Serra, 1,2 David Sanchez, 1 and Rosa Lopez 1 

1 Departament de Fisica, Universitat de les Illes Balears, E-07122 Palma de Mallorca, Spain 
2 Institut Mediterrani d'Estudis Avancats IMEDEA (CSIC-UIB), E-07122 Palma de Mallorca, Spain 

We analyze the spectral and transport properties of ballistic quasi one-dimensional systems in the 
presence of spin-orbit coupling and in-plane magnetic fields. Our results demonstrate that Rashba 
precession and intersubband coupling must be treated on equal footing for wavevectors near the 
magnetic field induced gaps. We find that intersubband coupling limits the occurrence of negative 
effective masses at the gap edges and modifies the linear conductance curves in the strong coupling 
limit. The effect of the magnetic field on the spin textured orientation of the wire magnetization is 
discussed. 
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I. INTRODUCTION 

Controllable manipulation of electron spins with elec- 
tric fields (gates) is a central requirement to spintronic 
devices. Semiconductor heterostructures offer the possi- 
bility of electric control of spins through intrinsic spin- 
orbit interactions. A major contribution to spin-orbit 
effects in two-dimensional (2D) electron gases of narrow- 
gap semiconductors originates from the macroscopic elec- 
tric field confining the electron gasi This implies an 
asymmetry in the quantum-well potential (the Rashba 
effect), which can be further tuned with gate voltages. 
The resulting spin-orbit coupling has been demonstrated 
in experimentsiSi&i 

In this paper, we are concerned with the Rashba 
interaction in ballistic quantum wires. Interestingly, 
these systems have been proposed as basic elements for 
practical applications such as spin-dependent field-effect 
transistors^ and spin filters^ and have been also consid- 
ered as detectors of entangled pairs of electrons^ and of 
the hyperfine nuclear spin dynamics. 8 The Rashba inter- 
action is described by the Hamiltonian, 



Oil 



-PxO-y 



a 2 



-PyO-x 



(1) 



where p — (p x ,Py) is the linear momentum and o~ x , o~ y 
are Pauli matrices. Spin transistors exploit the Rashba- 
induced precession of spins (a\ term for transport along 
the cc-direction) . However, in Eq. Q there is an addi- 
tional term, proportional to oli , which mixes nearest sub- 
bands with opposite spins and induces anticrossings at 
the degeneracy points of the wire's energy spectrumAifii 
Obviously, in a device one has a\ = a 2 , but it is of con- 
ceptual interest to distinguish in Eq. (JJJ between the two 
contributions (see below). 

When an in-plane magnetic field is externally ap- 
plied or it arises from ferromagnetic leads attached 
to the wire, the Zeeman splitting opens gaps in the 
wire spectrunij**2iiiiiaiiii strongly affecting the transport 
properties. It also influences the spin dynamics since 
whereas the Rashba precession randomizes the spin di- 
rection (a common spin quantization axis can be defined 
only for large wavevectors), the magnetic field tends to 



align the spins parallel to it. Many works negledi&iiii 2 ^ 
the effect of Rashba intersubband coupling (RIC) in 
quantum wires in the presence of in-plane magnetic fields. 
Here we demonstrate that both Rashba precession and 
intersubband coupling must be treated on equal footing 
for wavevectors near the subband gaps. We find that RIC 
hinders the formation of subband maxima, smoothing the 
gap edges and, as a consequence, strongly affecting the 
conductance steps. This conclusion is most relevant to 
the strong coupling limit, when spin-orbit and confine- 
ment energy scales are of the same order of magnitude. 
Section II presents the model and analyzes its charac- 
teristic energy subbands. In Sec. Ill we show the spin 
expectation values and magnetization distributions while 
Sec. IV focusses on the linear conductance and, finally, 
the conclusions are contained in Sec. V. 



II. HAMILTONIAN AND ENERGY SPECTRUM 

We consider a quantum wire, formed when a 2D gas 
is further confined in one direction [see Fig. ^ . The con- 
finement is assumed parabolic in the y-direction, Jt? C onf = 
mwoJ/ 2 /2, giving a wire orientation along x. An in-plane 
magnetic field B = (B cos 6, B sin (9) acts through the 
Zeeman Hamiltonian J^fz — gfJ-BB(cos9a x + sm9a y ) /2. 
Orbital magnetic effects are absent in this geometry since 
they arise only from perpendicular fields 4*^*1^ Adding 
all contributions, the resulting quasi-lD Hamiltonian 
reads, 



J?=(P 2 X + Pi)/"*™ + Jf conf + Jf z +Jf R . (2) 



Since the system is translationally invariant along x, 
the wave function dependence on this variable is e zkx , 
with k the propagation wavevector. It is then useful to 
recast the Hamiltonian for a given wavevector, J^, in 
dimcnsionlcss form^ 
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with the characteristic lengths Iq = ^h/mujQ (confine- 
ment), lz = y/h 2 /mg^sB (Zeeman coupling, g > 0)A& 
and U = h 2 /2mai (i G 1,2, Rashba interaction). Wc 
use in Eq. J2J the bosonic operator a k (at.) which lowers 
(raises) a subband index for fixed k, thereby n k = a k a k is 
the number operator. The third term in the right-hand 
side of Eq. © describes the free propagation motion in 
the x-direction whereas the forth and fifth contributions 
correspond to Rashba precession (h) and to RIC (h), 
respectively. The h Rashba coupling induces subband 
spin splitting whereas it is clear from Eq. that the 
I2 term couples adjacent subbands with opposite spins. 
In writing Jff k we have omitted the Coulomb interac- 
tion between the electrons since its effect can be taken 
into account, at least in part, via a renormalized Rashba 
coefficient iii 

Due to the presence of RIC, J%% is not diagonalizable 
except for a special wavevector (see below). Yet, it is easy 
to show that for fields parallel to the wire the two Rashba 
terms contribute equally. Without Rashba interaction 
the energy dispersion of a given subband is = (n + 
1/2) + (kl a ) 2 /2 + i 1 (l /l z ) 2 /2 with n = 0,1,2,... and 
7] = ±1. We find the energy correction to first order in 
the Rashba couplings: 
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Equation (0J already shows two noticeable features: (i) 
When 9 ^ there is a Rashba-induced splitting propor- 
tional to fc, well known from studies of spin-orbit effects 
at B — (see, e.g., Ref. fish . which combines with the 
Zeeman splitting to yield multiple subband crossings; and 
(ii) RIC (the h term) does not contribute to first order 
but for 8 — the Rashba precession (h) term is also 
zero. This implies that when 9 = both Rashba terms 
contribute to second order in the couplings. The full 
second-order correction reads: 
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FIG. 1: (Color online) Schematic representation of the wire 
potential Jtf? CO n{ = V(y) = mujQy 2 /2 and magnetic field orien- 
tation considered in this work. 




FIG. 2: (Color online) (a) Dispersion relation for = 0, lz = 
2lo, and h = Iq. Solid line corresponds to the case where the 
Rashba intersubband coupling is fully included (I2 — h) while 
dashed line shows the case without it (h = 00). (b) Same as 
(a) for = 7r/2. 



Adding all corrections one has En^ — E, 

p(2) 
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^nfcTj' wmcn is valid for l ,l z < h,h and l ^ l z . It 
should be emphasized that in the above perturbative 
analysis we have taken as expansion parameters lo/h 
and lo/h, assuming that the remaining factors of the two 
Rashba terms are similar, i.e., that klo — 1- Obviously, 
in the limit of very large k the Rashba precession term, 
proportional to klo, will clearly dominate over the RIC. 

Further progress is made if we diagonalize the total 
Hamiltonian in the absence of RIC, i.e., setting h — > 00 
in Eq. The exact eigenstates and eigenenergies (in 
hujo units), 
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are labeled, as before, with three quantum numbers: the 
propagation wavevector k 7 the index n of the correspond- 
ing ID harmonic oscillator wavefunction 4> n (y), and the 
branch-splitting quantum index rj = ±1. The nota- 
tion has been simplified by defining the complex num- 
ber z k = (lo/lzfeT 16 /2 + i(kl /2){l /h) and its argu- 
ment ri/j = Arg[zfc]. This complex quantity is taking into 
account the relative importance of Zeeman and Rashba 
precession energies, as can be seen more clearly when 
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rewriting it like z k — [(g/isB /2)e~ lS + ia\k\/hu>o. 

The energy spectrum of Eq. J7J), for 9 = and 9 = tt/2 
is plotted with dashed lines in Figs. E{a) and (b) for 
lz = 2Zo and strong spin-orbit l\ = Iq (taking a rj 10 
meVnm and InAs parameters the above values would 
correspond to Hluq rj 0.2 mcV and B w 0.1 T, which 
can be achieved experimentally)^ For both angles the 
spectrum shows energy crossings between subbands with 
n 7^ n! and opposite 77. For 9 = tt/2 there are addi- 
tional crossing points for klo ~ 0.2 between subbands of 
opposite r\ and equal n. Noticeably, for 9 = the Zee- 
man field produces gaps at k = and around that point 
the spectrum shows local maxima for the rj = — 1 sub- 
bands. From Eq. Q we see that a local maximum (min- 
imum) occurs at the 77 = — 1 subbands when lz > \p2l\ 
(lz < V^h). Thus, changing the magnetic field affects 
dramatically the spectrum and, as discussed below, the 
transport properties. 

Solid lines in Figs.|2fa) and (b) show the effect of RIC. 
We have obtained the spectrum with a direct numerical 
diagonalization of Eq. © in the basis of Eq. ©, trun- 
cating to a large enough n. The intersubband matrix 
elements of J^ a2 = (iZo/2v / 2?2)(a^ — a)a x read 



(nkr]\J^ a2 \n'k'ri') = 



4^2 h 



x (S n ', n -iVn - 8 n ',n+Wn + 1 ) . (8) 

As an alternative method, we have also discretized 
in real space with finite differences and diagonalized the 
resulting matrix. Both methods give identical results. 
When RIC is fully included we find that the crossings of 
different-n subbands (dashed lines) become anticrossings 
(solid lines) as expected. Importantly, for 9 = and high 
energy subbands RIC converts the maxima into minima 
[see, e.g., the second and third subbands in Fig^a)]. 

Figure [21 also shows a conspicuous downward shift of 
the solid lines with respect to the dashed ones for k- values 
close to the central gaps (for 9 = 0) or degeneracy points 
(for 9 = 7r/2). This can be easily explained by noting 
that the minimum gap condition for 9 ^ tt/2 or crossing 
point for 9 = tt/2 is given by Im[zfc s ] = with solution 
kg = (h/lz) sin 6*. When this occurs the Hamiltonian © 
can be exactly solved since J^ a2 does not couple states 
with different 77-indices [see Eq. (JHJ) and note that 77+77' = 
277<5 W ']. The Hamiltonian for each 77, 
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represents a harmonic oscillator with a shifted transverse 
momentum and a global negative energy constant. As a 
result, the exact states at the minimum gap are those 
of Eq. (JHJ shifted in momentum, e - %r im.a 2 y/h ^ nk ^ ri ( x ^y^ 
and the exact eigenenergies show a rigid shift E n k v — 

il/sqii 



For 9 = tt/2 [Fig.J^b)] we find new energy degeneracies 
entirely due to the presence of RIC in addition to the 
anticrossings similar to the 9 = case. For < 9 < tt/2 
there is a local minimum or maximum at k g depending 
on the size of ai with respect to ct\ and B. Our results, 
thus, demonstrate that including intersubband coupling 
is essential in any theory of strong Rashba interaction in 
quantum wires in the presence of magnetic fields and for 
k values near the subband gaps. 

To see how sensitive the local maxima are to the ef- 
fect of RIC we calculate the second-order correction in 
Ji? a2 to the energy spectrum at 9 = and then find the 
dimensionlcss effective mass at the extremum of the n?7 
branch, m/m* nil = d 2 E nkll /d{kl ) 2 \ k=kg . The sign of m* 
determines the character of the extremum, electron-like 
(to* > 0) or hole- like (to* < 0): 



1 + ^ ^ 
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For a 2 = and large a% (h < lz/V2), to*^ becomes 
negative for 77 = — 1 subbands, leading to a maximum 
in the energy spectrum, in agreement with our previous 
conclusions. However, for strong enough RIC the effec- 
tive mass sign is reversed. We emphasize that this is an 
effect purely due to the RIC and it is even stronger for in- 
creasing n, confirming the numerical results of Fig. [2TaV 
Equation (|10f> is a perturbative result and works well 
provided the confinement is the smallest length scale 
(la < lz, h,h)- Only for very wide wires (h,h, h < la) 
we numerically recover maxima (to*^ < 0). 



III. SPIN ORIENTATION AND 
MAGNETIZATION 

In general, when both Rashba terms are present spin 
is not a good quantum number and we find spin tex- 
tures, with the spin direction depending on k and the 
wire transversal coordinate y. On the contrary, if RIC 
is neglected the states are proper eigenspinors even in 
the presence of an in-plane magnetic field^i as seen from 
Eq. JHJ. Therefore, the existence of clear spin textures 
is a signature of the RIC term. The local spin compo- 
nents for the lowest subband at three different propaga- 
tion momenta are shown in the upper plots of Fig. [3] 
While all in-plane spins are essentially collinear, a size- 
able z-component, similar in magnitude to the in-plane 
one, precludes the definition of a common spin axis when 
klo 7^ and thus shows the importance of RIC. The local 
z-magnetization in real space (<J z (y)) is antisymmetric in 
y, leading to a vanishing integrated (a z )^ and giving 
rise to spin accumulations at the wire edges which are 
reminiscent of the intrinsic spin Hall effect but here 
the effect arises in a confined system^. 
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FIG. 3: (Color online) Lower plot: Dependence with k of the 
spin expectation values in the lowest subband. Solid lines are 
obtained including both Rashba terms while dashed lines are 
the results when RIC is neglected. We set lz = 2.2£o, l\ = 
I2 = lo and 9 = 0. Upper panels display the spin texture for 
three selected propagation momenta and taking into account 
both Rashba terms. The klo values (indicated by the thick 
arrows pointing on the fc-axis) are -0.75, and 0.75 for left, 
center and right upper plots, respectively. The vector plot 
shows the in-plane spin and the continuous line corresponds 
to the z component. 



There are exceptional cases where the spin is well de- 
fined even in the presence of both magnetic field and RIC. 
For |fcZo| S> 1 one always finds asymptotic eigenstates of 
cr y! 3I irrespectively of the magnetic field orientation 6, 
as seen in the lower plot of Fig. [3] For states with large 
negative (positive) k the spin is (cr y ) = — 1(+1). Impor- 
tantly, for k — (the minimum gap point for 9 = 0) we 
numerically find (a x ) — 77, regardless of the strengths of 
magnetic fields and intersubband coupling. This is clear 
from Eq. © , which shows that the spinors arc eigenstates 
of a x at the minimum gap momentum. Therefore, there 
is always a given wavector k g — (li/l z )sva.8, satisfying 
f2/c =0, for which the propagation direction is also the 
spin quantization axis. Away from this special point, a x 
tends to zero (see Fig. [3J| more slowly for stronger fields. 
Lower plot in Fig. [3] displays also the spin expectation 
values (o~ x ) and (a y ) when RIC is not taken into account 
(dashed lines). As we will see later, the effect of RIC in 
higher subbands is more notable but we emphasize that 
even in the lowest subband the effect of Rashba intersub- 
band coupling is not negligible due to the observation of 
spin textures. 

We plot in Fig. 0] the spatially dependent magnetiza- 
tion and the spin expectation values for the second sub- 



FIG. 4: (Color online) Same as Fig. 3 for the second subband. 



band with the same parameters as in Fig. [3] The effect 
of RIC is more pronounced than in the lowest subband as 
seen from the spin textures, which contain a noncollinear 
distribution even for the in-plane spin. In addition the 
local z-magnctization is bigger than in the lower sub- 
band and exceeds the horizontal components, leading to 
large spin accummulations at the wire edges. For k = 
we again find that all local spins are properly aligned in 
the direction of transport. Remarkably, the spin expec- 
tation values (a x ) and (a y ) change drastically when the 
Rashba interaction is fully included as shown in lower 
Fig. When RIC is absent, both expectation values dis- 
play abrupt jumps due to the subband crossings. As we 
pointed out earlier, RIC avoids these crossings and this is 
reflected in a much smoother behavior of (a x ) and (cry). 

Opposite sign accumulations of the z-magnetization at 
the wire edges are already found in the B = casei 19 i 20 
Differences with respect to the zero field case, however, 
can be summarized in that the magnetic field causes 
nonzero x-magnetization across the wire for wavevectors 
around the minimum gap point. The tilted distributions 
of in-plane magnetizations in Fig. 3 are due entirely to 
the ^-component of the field. Additionally, the texturing 
of the horizontal magnetization of the second subband 
(Fig. 4) is due to the combination of B and RIC. The 
magnetic field also shifts the minimum gap momentum 
kg from its zero value at vanishing B, for which an exact 
solution was obtained in the preceding section. 
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FIG. 5: Variation of the conductance with Rashba intensity 
and the Fermi energy fi in the full Rashba model (h = h, 
upper panel) and neglecting RIC (fo = oo, lower panel). We 
set lz = 1.8Zo and parallel magnetic field 9 = 0. A rigid shift 
with Rashba intensity has been taken into account defining 
(1/Hujo = u/fiujo + lo/Sli- The numbers in the plateaus give 
the conductance in units of the conductance quantum. Right 
small plots show vertical cuts of the corresponding left figures 
for the given values of Iq/21\. 



IV. TRANSPORT PROPERTIES 

The linear conductance Q at zero temperature is di- 
rectly related to the wire subband structure through the 
formulai2iii 



is generally reduced when both terms are included. 

For 9 = we plot in Fig. [S] the linear conductance for 
the cases with RIC, 1% — l\ (upper panel), and without 
it, l<z = oo (lower panel). When RIC is not included, 
Q alternates, with increasing fi, steps of +2e 2 /h with 
downward jumps of —e 2 /h due to the presence of max- 
ima in the lower subbands When RIC is included two 
major modifications are apparent: a) Stronger Rashba 
couplings are needed to observe the alternate steps of 
+2e 2 /h and —e 2 /h, and b) in any case this anomalous 
pattern of steps vanishes when increasing the Fermi en- 
ergy. Particularly, as shown by the insets, when Zo/2/i = 
1 the sequence of conductance steps in units of Qq is 
+2, — 1, +1, +2, — 1, +1, . . . when RIC is neglected, and 
+2, — 1, +1, +1, +1, . . . when it is included. Therefore, 
Fig. Improves that to observe the modifications of the lin- 
ear conductance steps due to the Rashba interaction it is 
essential to have a relatively low Fermi energy or, equiv- 
alently, a rather small number of propagating modes. 

Figure shows the case 9 = ir/2. In lower Fig. El when 
I2 = 00, the conductance consists of regular jumps of 
one conductance quantum. From Fig- 13b) we determine 
that the width of the odd conductance plateaus corre- 
sponds to the Zeeman energy whereas the width of the 
even plateaus amounts to 1 — (lo/lz) 2 (in units of huio). 
This physical scenario is strongly modified when RIC is 
taken into account, as seen in the upper Fig. where the 
conductance shows a much richer behavior. It is remark- 
able the occurrence of restoration of the spin degeneracy 
for special values of the Rashba coupling and magnetic 
field where the third and fifth conductance plateaus col- 
lapse. In general, the tendency is to suppress regions for 
odd values of Q . Thus, Q presents wide (narrow) plateaus 
for even (odd) values of the conductance quantum but 
their widths are not trivially related to either the Zee- 
man energy or the confinement energy alone. Besides, at 
large Rashba strengths we also find decreasing jumps in 
the conductance. 



^=TE e ^-^ ) )sgn(mp) 



(11) 



CONCLUSIONS 



which accounts for the number of occupied subbands 
propagating in a given direction and the conductance 
quantum per spin component, Qq = e 2 /h. In Eq. /i 

is the Fermi energy while £^ and mf'* denote, respec- 
tively, the energy and effective mass at the j-th local ex- 
ternum in the i-th subband. We have also defined <d(x) 
and sgn(x) as the step and the sign functions. For sub- 
bands containing a single minimum Q only attains e 2 /h 
steps when [i is increased. In the case of more compli- 
cated subbands it is easy to see from Eq. that the 
linear conductance step is doubled whenever the Fermi 
energy exceeds two local minima but is below the en- 
ergy of the maximum lying in between them. Once fj, 
exceeds the energy of the maximum, the conductance is 
decreased by e 2 //i&ii As discussed above the RIC can 
convert subband maxima into minima and, therefore, Q 



We have discussed the spectrum, spin orientation and 
linear transport of ballistic quantum wires in the pres- 
ence of Rashba interaction and in-plane magnetic fields 
at an angle 9 with the propagation direction. Analytical 
and numerical results demonstrating the crucial impor- 
tance of the Rashba intersubband coupling for all these 
properties and for k values near the subband gaps have 
been presented. At strong spin-orbit coupling the modi- 
fications in the energy subbands include the appearance 
of sizeable shifts, anticrossings as well as big reductions 
of the subband maxima. The changes in subband struc- 
ture lead to precise predictions for measurements of the 
linear transport properties, such as a severe reduction of 
the anomalous conductance steps when 9 = and a non- 
trivial dependence of the steps with the Rashba intensity 
for 9 = 7r/2, with collapsing points for the odd plateaus. 



6 




Spin textures in the local magnetization are obtained 
only when the Rashba intersubband coupling is taken 
into account. Opposite-sign accumulations of the spin z- 
magnetization at the wire edges are strongly dependent 
on the fc-value and subband index. For subbands above 
the lower one we have shown that the horizontal magne- 
tization is also textured. At a particular momentum (the 
minimum-gap momentum), depending on magnetic field 
strength, orientation, and Rashba intensity, the analyti- 
cal solution having the spin along the transport direction 
has been obtained. We believe these results are relevant 
for the future design of spin transistors and spin filters. 



FIG. 6: Same as Fig. 5 for perpendicular field (9 — 7r/2). 
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